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Abstract 

We consider the indeterminacy in the sign of the neutrino asymmetry generated by active- 
sterile neutrino oscillations in the early universe. The dynamics of asymmetry growth is 
discussed in detail and the indeterminacy in the final sign of the asymmetry is shown to be a 
real physical phenomenon. Recently published contradicting results are carefully considered 
and the underlying assumptions leading to the disagreement are resolved. 



The recent observation of strong zenith angle dependence of the atmospheric neutrino 
deficit by the Super-Kamiokande neutrino experiment has provided strong evidence for neu- 
trino oscillations — ux, where i>x is either Uj. or a new, sterile neutrino |l|]. While the 
^'/x ~ i^s solution presently is less favored by SK data [Q, reconciling the existing data from 
all the neutrino experiments, including LSND, is not possible unless there exists at least 
one sterile neutrino mixing with the active neutrinos. Such mixing would have interesting 
consequences for primordial nucleosynthesis and CMB radiation |T^. For example, 

sterile neutrinos could be brought into equilibrium prior to nucleosynthesis, increasing the 
energy density of the universe and thereby neutron-to-proton freeze out temperature, leading 
to more helium-4 being produced. This scenario has been numerically studied and strong 
limits to neutrino mixing parameters have been obtained 

Under certain conditions active-sterile neutrino oscillations may also lead to exponential 
growth of neutrino asymmetry, as was first discovered by Barbieri and Dolgov P] . Later Foot 
and Volkas observed that by this mechanism very large asymmetries could be generated, 
which would have a significant effect on the primordial nucleosynthesis by directly 

modifying directly the n ^ p reactions. Moreover, they showed that an asymmetry generated 
by Ut- — Ug mixing could suppress sterile neutrino production in i/^ — z/^-sector, loosening the 
bounds of 0], according to which the SK atmospheric deficit could not be explained by 
z/^ - Ug mixing. 

Later it was found that this asymmetry generation is chaotic in the sense that determining 
the sign of the final asymmetry sign{L) does not simply follow from the initial conditions [|lT 



This phenomenon was studied in [|T^ , and it was shown that the indeterminacy is associated 
with a region of mixing parameters, where asymmetry is rapidly oscillating right after the 
resonance. As a consequence the the amount of Helium-4 produced cannot be precisely 
determined in such a scenario 1112 . 



In a recent paper |^ it was claimed, however, that sign{L) is completely determined 
by the initial asymmetry, and moreover, that there is only a slight growth of asymmetry 
after the resonance. In this article we clarify the origin of indeterminacy in sign{L) and 
show that it is a real physical phenomenon, not disturbed by numerical inaccuracy. Instead, 
we will argue that the disagreement arises due to overly simplifying approximations used 
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in [^. Finally, we will point out to a likely cause leading to observed suppression of the 
final magnitude of the asymmetry in ||13|| . 



In the early universe the coherent evolution of the neutrino states is interrupted by 
frequent decohering collisions. Therefore the evolution of the system needs to be studied 
using the density matrix formalism. We parameterize the reduced density matrices of the 
neutrino and anti-neutrino ensembles as 

p,^^Po(l + P), Pp = ^Po(1 + P), (1) 

where each matrix is assumed to be diagonal in momentum space, while each momentum 
state has 2 x 2-mixing matrix structure in the flavour space. Solving the full momentum 
dependent kinetic equations for these density matrices numerically n] is a very compli- 



cated task and all attempts published to date have used some approximations to simplify 
the problem. Here we use momentum averaged approximation i.e. we set P{p) P {{p)), 
with (p) ~ 3.15T. This approach has been found to give a very good approximation of 
the equilibration 0], and it will be sufficient for the purposes of this letter. The coupled 
equations are then (in the case of z/^ — oscillations, other cases can be obtained easily by 
simple redefinitions which are found for example in |^): 

p = VxP-(Z} + ^logPo)PT+(l-P.)^logPo2, 
P = VxP-(D + ^logPo)PT+(l-n)^logPo2, 



Po = (r(z/^z/r (nig - riy^riy^j , 

Po = (T{yr^T aa)) ("-eg - n^^n^r) , (2) 

where x = dx/dt and P^ = PrX + PyY- The damping coefficients for particles and anti- 
particles are ~ Z) ~ 1.8G|^T^ very accurately. The rotation vector V is 

V = K-x+(H + V^l)z, (3) 
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where 

V. = ^sin2^, 
2{p) 

Vo = —-cos2e + 5Vr, 
2{p) 

Vl = -V2GfN^L, (4) 

where 9 is the vacuum mixing angle, = ml^ — ml^, is the photon number density 
and the effective asymmetry L in the potential Vl is given by 

L = -^Ln + L,^ + L,^ + 2L,XP) =ri + 2L,^{P), (5) 

in the case of an electrically neutral plasma. Asymmetry L^^ is obtained from 

L,^ = I (Po(l + Pz) - Po(l + Pz)) (6) 



and Ln is the neutron asymmetry. The potential term 6Vr is approximately 

5K = n.SGpN,^^. (7) 

The rotation vector for anti-neutrinos is simply V(L) = V(— L). 

Neutrino and anti-neutrino ensembles are very strongly coupled in Eq. (Q) through the 
effective potential term Vl{L), which makes their numerical solution particularly difficult; 
as long as neutrino asymmetry remains small, there is a large cancellation in the Eq. (|^), 
leading to a potential loss of accuracy. To overcome this problem we define new variables 

P± = P^± P„. (8) 

In terms of these the Eq. (§) for becomes 

P± = -VoP^-VlP^-DP^, 
P± = VoP^ + VlP^-DP^^-V,P^, 

7± 



p^ = v;p,± + A± 2-P/ -A^P±, (9) 
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where we have defined 



d 
It 

A_ = 0, 



= -logPo+, 



b = D + A+. (10) 

Finally, since for the averaged interaction rates (F) = (f), the difference Pq is not affected 
by collisions and we find 



Po = 0. (11) 



Our objective is to study whether the sign of the final asymmetry, sign{L) , follows 
deterministically from the initial conditions. In |TT],0| it was found to be chaotic, whereas 



the authors of []T^ claim that sign{L) is fully deterministic and equal to the sign of the initial 
neutrino asymmetry. Similar results have been reported by other groups as well f^^j • 

The key ingredient in the physics leading to the growth of the asymmetry is the appear- 
ance of the resonance. Indeed, if the squared mass difference 6m'^ < 0, the effective potential 
Vo±Vl goes through zero at the resonance temperature 

Te ~ 16.0 {\6m^\ cos 2^)^/^ MeV, (12) 

where ~ is assumed, as effective asymmetry L is driven to zero well before the resonance. 
After the resonance the balance of the system abruptly changes which leads to a rapidly 
growing L. We have numerically solved Eqs. (^ and ([TT|), and examples of the results are 
shown in Fig. |I|. 

The behaviour of the system can be understood by a simple analogy of a ball rolling 
down a valley. After the resonance temperature Tc the originally stable valley at L = 
becomes a ridge line separating two new, degenerate valleys corresponding to solutions of 
Vo ± Vl = 0. The system may first oscillate from one valley to another passing over the 
ridge line. However, because of friction (represented by damping terms) it will eventually 
settle into one or the other of the new valleys. It is easy to imagine that when there are 
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many oscillations, even a very small difference in initial conditions may grow to a large 
phase difference at the time of settling down. In Fig. Q these effects are demonstrated for 
two initial values: t] = 10~^° (solid line) and 77 = 2 x 10^^° (dashed line) and oscillation 
parameters 5m^ = — 10~^ eV^ and sin^ 26' = 10~^'^. For these parameters the resonance 
occurs at temperature = 7.41 MeV. 
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Figure 1: The evolution of neutrino asymmetry at the resonance temperature for mixing 
parameters = — 10~^ eV^ and sin^ 2^ = 10~^'^. Solid line corresponds to initial value 
77 = 10~^° and dashed line to rj = 2 x 10~^°. 

Let us point out that change of sign{L) can be effected either by variations in the os- 
cillation parameters (5m^ and sin^ 29 or variations in the initial asymmetry t] as was shown 
in [|12|. These two cases are physically quite different of course. Varying oscillation param- 
eters changes the shape of valleys forming after the resonance. This is an important issue 
because there will always be some experimental uncertainty in the measurements of masses 



and mixings, leading to unavoidable uncertainty of SBBN predictions in this scenario [|T2 
Variations in rj are physical, for example due to local inhomogenieties in the baryon asymme- 
try created by early phase transitions, and correspond to deviations in the initial conditions 
(speed and direction of the ball) upon entering the resonance region. 



In a recent article by Dolgov et.al. |13|, results were presented which are in sharp dis- 



agreement with ours. In particular, it was claimed that the sign of neutrino asymmetry is 



completely deterministic. Moreover, they suggested that the chaotic behavior seen in []TT],|12 
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would be due to accumulated errors in the numerical codes. This suggestion is not correct. 
The numerical errors are completely under control in our computations. We have checked 
that allowing for the local error tolerance magnitudes larger than actually employed to get 
our results, have no effects on them. Rather we will now argue that the disagreement is due 



to ill-justified analytic approximations used in [|T3| to simplify the problem, which lead to 
artificial sign determinacy. 



First approximation made in was to neglect the term VxP~ in the equation for P~ 
under the assumption that it is small in comparision to V^Py^'-term. However, just before 
the resonance effective asymmetry L is driven to zero, so that effective potential Vl is in fact 
very small, and VxP~ should not be expected to be subdominant. We have indeed found 
that this term is of crucial importance for the initial asymmetry growth, as it prevents L 
from getting arbitrarily small value before the resonance, as will be discussed below. 



Second approximation made in has even more dramatic effects. There it was argued 
that because Eq. (P) for P~ and P~ can be scaled to a form P~y = Q{a + 6 + c), where 
Q ~ 5.6 X lO'^y^l cos266m'^\ is a large parameter, it should be safe to set the derivatives P~y 
to zero. The rapid oscillations seen in our solutions indicate this approximation breaks down 
at the resonance temperature. 

To study the effect of these approximations quantitatively, we introduce them into our 
equations. Dropping the term V^P' and setting P^^ to zero in the equations (§) we then 
find the constraints 

= -V0Py--VLPy^-DP-, 

= VoP- + Vi^P;^-DPy-. (13) 
From these equations, one can solve the evolution of P~ and Py algebraically with the result 

The remaining variables in Eqs. (^ and ( |TI| ) are then solved numerically. In Fig. ^ we plot 
the results of a computation with (solid line) and without (dashed line) the implementation 
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of the constraints (|T4|) for same parameters as in Fig. [^. The constrained solutions which fall 
on top of each others in the figure are indeed fully deterministic and display no oscillation. 
This is to be expected because the solutions ([I^) cannot produce sign changing oscillations 
in L, since when L goes to zero so does P~ and hence P~ is then strongly suppressed in 
the constrained case. This prohibits L from changing sign, so that sign{L) is decided by 
the initial value of L. When the evolution of the off-diagonal components is neglected, the 
mechanism of the asymmetry growth is different, and instead of oscillating between the two 
new valleys, the system slowly rolls down to one of them after the resonance. This is in 



fact common characteristic to so called static approximations [§|,|T3[! where the differences 
between off-diagonals in the neutrino and antineutrino density matrices are ignored. 
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Figure 2: Comparision of constrained (solid line) and unconstrained (dotted line) solutions 
for the same parameters as in Fig. |I] 

Let us now consider in detail how the complicated interplay between the various terms in 
the Eq. (^ leads to the initial exponential growth and oscillations of the asymmetry at the 
resonance temperature. Before the resonance P"*" components are practically independent of 
components as effective asymmetry L is very small. However, off-diagonal components 
and Py do grow to fairly large values. When Vq changes sign at the resonance P^ 
is rapidly driven to change sign, while the evolution of P^ remains unaffected due to the 
additional term V^P^. P^ stays near its initial value until the off-diagonals begin to grow, 
after which it begins to diminish, signalling the sterile neutrino production. Overall the 
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evolution of P"*" components is smooth (see Fig. 
resonance is the changing of . 



and the noticeable direct effect of the 
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Figure 3: Evolution of (solid line) and P^*" (dashed line) at the resonance. 



The evolution of P~ components is more complicated. Before the resonance neutrino 
and anti-neutrino ensembles follow each other closely in the sense that off-diagonals P~y 
are small. The main force which is keeping the difference of neutrino and antineutrino off- 
diagonals small is the potential Vl and not the effect of damping terms. This can be confirmed 
by explicitly plotting individual terms appearing in the derivatives P~y (see Fig. ^). The 
underlying mechanism driving P~y close to zero before the resonance is then the cancellations 
between the effective potential terms. The Kr-PT term in Eq. (^) prevents L from stabilizing 
to an arbitrarily small value, since when L and off-diagonal components P~y are driven 
towards zero before the resonance, the small difference P~ becomes important. The value 
of Vl together with differences P~y will cancel the effect of a non-zero P~ . In this way the 
terms VqP" and ViPy cancel each other in the equation for P~, and VqP' together with 
VlP^ cancel V^PT in the equation for P~. 

After Vo changes sign the cancellation of the potentials in P~ equation will not work, 
which turns around the effect of Vl in the P~ equation. Moreover, as P^ changes sign due 
to the resonance, the effect of Vl changes also in the P~ equation. The difference between 
both off-diagonal components is now growing due to a non-zero value of L leading to a self 
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Figure 4: The values of invidual terms in the Eq. (^) for Py during the resonance: Py (thick 
sohd hne), VlP^ (sohd hne), VqP" (long dashed line), VxP^ (short dashed line) and -DP" 
(dotted line) 



supporting exponential growth of L. We wish to stress that the magnitudes of all the effects 
discussed here are well above the numerical accuracy. It is true that before the resonance L 
goes very close to zero, but this is not relevant. What is relevant instead is that P~ remains 
at a value of order rj. 

It is interesting to note that the dynamics explained above will always lead to an initial 
growth of L into the direction given by the sign of rj. This results agrees with the analytical 
considerations in 0. Nevertheless, as we have seen, this does not guarantee that the final 
sign is that of rj. One should bear in mind that the beginning of the asymmetry growth 
is a very complicated phenomenon, where all the variables and almost all the terms in the 
evolution equations are important to the outcome of the resonance and a complete account 
of the dynamics of the off-diagonal elements is of crucial importance. This makes it of course 
very difficult, if not impossible, to find any sensible analytical approximation to Eq. (|^). 

There still remains one contradiction with respect to the results of |T^ and the constrained 
solution obtained here; namely, as seen from Fig. 0, even the constrained solutions yield L 
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which grows to a large value. The quantitative study of this effect is beyond the validity 
of the method employed here. We may however point out a further approximation made 
in 



13[| , which appears to be the reason why they see a much weaker asymmetry growth. The 



collision term in [|l^ for the active- active component, paa of the density matrix (0) is of the 
form 

r(p)(Paa(p)-/eg), (15) 

where in [0] f^q was taken to be the free Fermi distribution function f^q = [1 + exp(p/T)]^^. 
The authors in note themselves that feq should actually be the distribution function 
which includes a chemical potential, feqifJ^) = [1 + exp(p/T =f fi/T)]^^, but they assume the 
difference to be minor because the /i/T is small. However, one can show that fi/T ^ 0.7 L^^ 
and expand feq{L^^) to give 

UU,)«U0)±0.7L..^^M^^, ,16) 

Approximating the second term of the expansion to zero in ( ]TBD then means that the system 
is seeking the free Fermi distribution instead of the correct equilibrium distribution which 
includes an asymmetry. This gives rise to an artificial force proportional to r{p)L^^ resisting 
the growth of the asymmetry. Taking this into account, together with their previous as- 
sumptions leading to constraints (p!^ , according to which the rate of asymmetry generation 
itself is weak 

P- oc P- ocVlocL, (17) 

it appears likely that this force is able stop the asymmetry growth before the system has 
reached the true bottom of the valley. We chose not to pursue this issue further, since even 
the previous approximations has rendered the system unphysical by denying the possibility 
of chaoticity. 



In this letter we have considered the indeterminacy or chaoticity in the sign of the neutrino 
asymmetry L^^ arising from z/^ — Ug mixing in the early universe. We carefully discussed 
the dynamics giving rise to the growth of asymmetry, and unravelled the mechanism leading 
to the uncertaintity in sign{L) . We confirmed that in the region of the parameter space 
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identified in the system is very sensitive to small variations in the initial conditions 
and mixing parameters. We have carefully checked that our numerical methods are higly 
accurate, so that the effects of the variations leading to a sign indeterminacy in the final 
asymmetry are physical. 



As we pointed out, our results are in contradiction with the recent claims |T^, according 
to which the sign of asymmetry is completely deterministic and the asymmetry growth is 
small compared to the results of 10-^, |lT], |12| . We have showed that the results of |]13 



are 



in fact unphysical and that they arise because of oversimplyfying approximations which 
artificially stabilize the dynamics responsible for the chaoticity. 

Indeed, we found that all the terms in the evolution equations and all the components 
of the density matrix are important for the dynamics of the system. Hence it appears 
unlikely that any simplifying analytic approximation can be found that would describe the 
system adequately. We have employed the momentum averaged equations, as they are 
sufficient to study the chaoticity of the asymmetry growth and to resolve the validity of 
approximations imposed in [jl3l . Our unpublished results with a momentum dependent code 
support the results of this letter, as well as do the results obtained by other groups working 
with momentum dependent equations 0, |l8l, although the chaotic region appears to be 
somewhat smaller, when momentum dependence is included. 
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